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The growing probabilities of additional offspring with the beneficial reversal allele 
for various population sizes, N, sequence lengths, L, selective advantages, s, fitness 
parameters, k, and measuring parameters, C, were calculated in an asymmetric 
sharply-peaked landscape using the decoupled continuous-time mutation-selection 
model. The growing probability in the stochastic region was inversely proportional 
to the measuring parameter when C < 1/Ns*, bent when C ~ 1/Ns* and became 
saturated when C > 1/Ns*, where s* is the effective selective advantage. The 
saturated growing probability in the stochastic region was approximately the effective 
selective advantage when C 1/Ns* and s* <C 1. The present study suggests 
that the growing probability in the stochastic region in an asymmetric sharply- 
peaked landscape in the decoupled continuous-time mutation-selection model can be 
described using the theoretical formula for the growing probability in the Moran two- 
allele model. The selective advantage ratio, which represents the ratio of the effective 
selective advantage to the selective advantage, does not depend on the population 
size, selective advantage, measuring parameter, and fitness parameter; instead the 
selective advantage ratio decreases with increasing the sequence length. 
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I. INTRODUCTION 


Two mechanisms, random mutation and natural selection, are considered in the mutation- 
selection model to explain an evolutionary process. The mutation-selection model is gen¬ 
erally classihed into two types, the coupled model and the decoupled model, depending on 
the mutation mechanism considered. The coupled model considers mutations that are only 
allowed as replication errors during reproduction events. The decoupled model considers 
mutations that are allowed at any time during the life cycle, such as those due to the effects 
of radiation, free radicals or thermal fluctuations[l-5|. Hence, mutation and selection are 
assumed to be independent processes that occur in parallel. 

From a biological point of view, the relative contributions of the underlying mutation 
mechanisms for the two versions of the mutation-selection model are controversial. A muta¬ 
tion mechanism described by the coupled version implies a relatively constant mutation rate 
per generation, whereas the mutation mechanism described by the decoupled version implies 
a .datively coaetaat .atatioa .ate ia t....e. Cow aad Ki...u.afl and Akinj^ studied the 
decoupled mutation-selection model. Baake et al. reported that the decoupled mutation- 
selection model was equivalent to an Ising quantum spin and calculated the observable 
quantities in three htness landscapes using the methods from statistical mechanics js-S]. 
Recently, Saakian et a/.jbl and Park and Deem[^ examined the decoupled continuous-time 
mutation-selection model using a functional integral representation. 

Several studies have considered the evolutionary process transiting to a neighboring higher 
htness peak through the low-htness valley. Iwasa et al. calculated the rate of stochastic 
tunneling in a large population when the intermediate mutant had a lower htness than the 


wild type, but the second mutant was benehcial[8|. Weissman et al. provided a complete 
description of the rate at which an asexual population traversed a specihc path through 
a genotype space to a particular htter genotype by neglecting backward mutations and 
genetic diversity}^. Biological examples of the low-htness valley can be found in pathogens 


that escape the immune response of their hosts, such as Escherichia coli 
typhimurium 12f| and the human immunodehciency virus 




Salmonella 


, Another example might be 

the types of cancer that occur only after a series of mutations[l^. 

Previous studies calculated the crossing time and the growing probability in the cou¬ 
pled discrete-time mutation-selection model by switching on an asymmetric sharply-peaked 
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landscape, such that the reversal allele of the optimal allele had higher htness than the op 


timal allele, from an initial steady state in a sharply-peaked landscape 
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18| . The growing 


probability was dehned as the probability that additional offspring with the reversal allele 
grow to a measuring criterion. The present study calculated the growing probability in the 
decoupled continuous-time mutation-selection model, which was evaluated as the ratio of 
the hrst arrival time to the hnal arrival time. The hrst arrival time was dehned as the time 
when the hrst additional ohspring arrived at the reversal allele. The hnal arrival time was 
dehned as the arrival time of a successful ohspring with the reversal allele that could grow 
to a measuring criterion. 

The growing probability of additional ohspring with the reversal allele in the decoupled 
continuous-time mutation-selection model was calculated systematically using computer sim¬ 
ulations for various population sizes, sequence lengths, selective advantages, htness param¬ 
eters, and measuring parameters, and was compared with the theoretical formula for the 
growing probability in the Moran two-allele model 19|, l20|. This paper is organized as fol¬ 
lows: Section II briehy describes the decoupled continuous-time mutation-selection models 
and the computer simulation method. Section III presents the computer simulation results 
for the growing probability in the decoupled continuous-time mutation-selection models. 
Section IV summarizes the results. 


II. MODEL AND METHOD 

This study considered the decoupled continuous-time mutation-selection model in a self- 
reproducing haploid population of N individuals, whose size remains constant over the 
generations. The multiple alleles. Si, might be represented by a string of L spins, which use 
two values, -|-1 and —1. The time dependence of the number, rii, of individuals with the 
allele. Si, in the decoupled continuous-time mutation-selection model can be described by 
the following equation: 

hi{t) = {wi - w) ni{t) + ^ Qij nj{t) , (1) 

j 

where Wi is the fitness of the allele Si, and w{= is the mean htness of the 

population. In Eq. ([1]), the mutation matrix, Qij, which gives the mutation rate from Sj to 
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Si, excludes double mutationsjd, 5|: 


-L/i , 

if i= j , 



if dij = 1 , 

(2) 

0, 

otherwise , 



Qij ^ 


where fi denotes the mutation rate per sequence element per unit time, and dij denotes the 
Hamming distance between Si and Sj. 

Consider a sharply-peaked landscape such that the optimal allele Sq = (-1-1, -|-1, • • • , -|-1) 
has a maximum fitness 1, and the remaining alleles have a lower fitness, 1 — kL. The 
parameter k is called the htness parameter. Equation ([1]) was solved using a computer 
simulation with the assumption that all individuals had the optimal allele initially. The 
population then approached a stochastic steady state, which is a quasispecies near the 
optimal allele in a sharply-peaked landscape. After arriving at a stochastic steady state, the 
sharply-peaked landscape was changed suddenly to the following asymmetric sharply-peaked 
landscape: 


if s = s 


w = < 


0 ) 


^0 5 


( 3 ) 


1 , 

1 + s , if S = S, 

1 — kL , otherwise , 

where the maximum fitness was assigned to the reversal allele, Sq = (—1, ■ • • ,—1), and 
the fitness of the remaining alleles was unchanged. All other parameters were fixed in this 
procedure, and the parameter, s, in Eq. ([3]) is called the selective advantage of the reversal 
allele over the optimal allele. The reversal allele had higher fitness than the optimal allele, 
assuming that the selective advantage, s, is a positive value. The asymmetric sharply-peaked 
landscape had two fitness peaks that were separated by a low-fitness valley with a tunable 
depth and width. 

The numbers, n*, were grouped into L + 1 distinct classes according to their Hamming 
distances from the optimal allele, Sq, to improve the speed of the computer simulation 2l[ |. 
The number, Ni{= individuals in class I was defined by the sum of all numbers 

Hi with a Hamming distance, d = I, from the optimal allele. The kinetic equation, Eq. o. 
for the decoupled model can be expressed using the newly dehned Ni as follows: 


Niit) = {Wl - W) Niit) + ^ Tira Nmit) , 


m=0 


( 4 ) 
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Tim — < 


(5) 


—Lfi , if / = m , 

(L — m)fi , if Z = m + 1 , 
mfi , if Z = m — 1 , 

0 , otherwise , 

where Tim gives the mutation rate from class m to class Z. The relative density, Xi, of class 
Z was then calculated by dividing the number, Ni, by the population size. 

The computer simulations used the logistic branching process with the replacement of 


Ni{t) in Eq. (g]) by {Ni{t +dt) - Ni{t))/dt 


19 


22 


24| . The simulated population was evolved 


using the time steps of dt = 0.1 generations in this study. The mean fitness, w, of the 
population was calculated at the beginning of each time step. The branching process assumes 
that the population in the next time step is determined by the following two events: each 
individual in a class can mutate into any class with a certain probability, and each individual 
in a class can produce an offspring in the same class with a certain probability. More details 
of the simulation for a hnite population can be found in Refs. 19 and 24. The time was 
set to Z = 0 when the htness landscape changed suddenly to an asymmetric sharply-peaked 
landscape. The increase in the relative density AXd(t){= Xd(t) —Xd(t = 0)) was calculated. 
The computer simulation was completed when the increase in relative density with the 
reversal allele, AX^lt), achieved the measuring criterion AXL{t) = C. The parameter, C, 
is called the measuring parameter. 

The computer simulation recorded all arrival times, Z, satisfying two conditions, AXiit — 
1) < 0 and AXi(Z) > 0. The hrst arrival time, Zq, was dehned as the earliest among the 
recorded arrival times, and the hnal arrival time was dehned as the latest after the computer 
simulation is complete. Therefore, the hnal arrival time means the arrival time of a successful 
ohspring with the reversal allele that can grow to a size, NC. The growing probability of 
additional ohspring with the reversal allele was evaluated as the ratio of the hrst arrival time 
to the hnal arrival time, assuming that the additional ohspring with the reversal allele are 
supplied at average time intervals of Z^ by stochastic huctuations. This assumption holds 
well in the stochastic region, where the arrival of additional ohspring with the reversal allele 
is a rare event with a quite long interval. Fortunately, the main consequence of this study 
was obtained from the simulation results for the growing probability in the stochastic region. 
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III. RESULTS 


The computer simulation results for the growing probability were obtained by averaging 
M independent runs. The mean and the standard error were calculated from M data values 
transformed into a logarithmic scale because the data showed a Gaussian distribution on 
a logarithmic scale. The number of runs was set to M = 3000 throughout this study. 
The computer simulation results for an infinite population were calculated using the same 
method reported elsewhere 2^. Hereafter, AX£°(f) denotes the computer simulation result 
in an infinite population and AXi(t) denotes the computer simulation result in a finite 
population. Figure [H shows the logarithm of the growing probability, log(Pg), as a function 
of the logarithm of the measuring parameter, log(C'), for L = 10, fi = 0.001, k = 0.003 and 
s = 10“^^. In Fig. [U the simulation results for the growing probability in the decoupled 
continuous-time mutation-selection model are represented by the solid lines with symbols, 
which correspond, right to left, to the population sizes N = 10^, 10®, 10^, 10® and 10®. 

A previous study showed that the boundary between the deterministic and stochastic 


regions in the decoupled continuous-time mutation-selection model cou 


using the critical population size, Nc, satisfying NcAX'^{t = 0) = 1 2^. The computer 


d be determined 


simulation for an inhnite population showed that AX£°(f = 0) Ri 4.73 x 10“® and Nc ~ 
2.12 X 10® for L = 10, /i = 0.001 and k = 0.003. A previous study reported that the 
growing probability in the stochastic region in the coupled discrete-time mutation-selection 
model was inversely proportional to the measuring parameter when C < 1/Ns, bent when 
C Ri 1/As, and saturated when C > l/AsQ. 

Figure [H shows that the simulation results for the growing probability for N = 10®, 10® 
and 10^, which belong to the stochastic region because N <C Nc, are inversely proportional 
to the measuring parameter when C < 1/Ns*, bent when C ~ 1/As*, and saturated when 
C > 1/As*, where s* is the effective selective advantage. The simulation suggested that the 
effective selective advantage, s*, is a hundred times larger than the selective advantage, or 
s* Ri 100s Ri 10“^. Let the selective advantage ratio, r, be the ratio of the effective selective 
advantage to the selective advantage(r = s*/s). The simulation results suggested r r: 100 
for the parameter region in Fig. [1] The simulation indicated that the number of additional 
offspring with the reversal allele, ul, in the stochastic region drifts neutrally until reaches 
the size, NC, when ul < 1/s*, and that it grows deterministically when ul > 1/s*. Figure 
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[Ualso shows that the saturated growing probabilities, P*, for N = 10^, 10® and 10^ become 
approximately the effective selective advantage, s* ~ 10“^, when C ^ 1/iVs* and s* <C 1. 

Figure [2] shows the logarithm of the growing probability, log(Pg), as a function of the 
logarithm of the measuring parameter, log(C'), for L = 10, /i = 0.001 and k = 0.003. In 
Fig. O the simulation results for the growing probability in a hnite population with a size, 
N = 10^, in the decoupled continuous-time mutation-selection model are represented by the 
symbols, which correspond, top to bottom, to the selective advantages, s = 10“^, 10“^, 10“® 
and 10“®. The parameters in Fig. |2] belong to the stochastic region, N Nc, because the 
critical population size, Nc, is independent of the selective advantage. Figure [2] clearly shows 
that the growing probabilities for various selective advantages are inversely proportional to 
the measuring parameter when C < 1/Ns*, bend when C ^ 1/Ns* and become saturated 
when C > 1/Ns*, where s* ~ 100s. Figure [2] also shows that the saturated growing prob¬ 
abilities, P*, for various selective advantages become approximately the effective selective 
advantage when C S> 1/Ns* and s* <C 1. 

The following theoretical formula for the growing probability, Pg, represents the proba¬ 
bility for a benehcial allele to grow to a size, NC, in the Moran two-allele model without 
mutations, when the selective advantage between two alleles is assumed to be equal to the 
effective selective advantage, s* jl9l. l20l|: 

s* (1 -1- s*)^^-^ 


PI = 


(1 -f- s*)^^ — 1 

if s* > 1, s*NC > 1 , 
if s* < 1, s*NC > 1 , 
if s*NC < 1 . 


( 6 ) 


1+s* 
S* . 
1 


NC ’ 

In Fig. 121 the theoretical formula in Eq. ([6]) with N = 10^ and s* = 100s is represented by the 
short-dotted lines, which correspond, top to bottom, to the selective advantages s = 10“^, 
10“^, 10“® and 10“®. A comparison of the symbols and the short-dotted lines shows that 
the simulation results for the growing probability, Pg, can be htted well using the theoretical 
formula for the growing probability, Pg, in the Moran two-allele model. A previous study 
showed that the growing probability in the stochastic region in the coupled discrete-time 
mutation-selection model could be described using the theoretical formula for the growing 
probability in the Wright-Fisher two-allele model[18[. The present study suggests that the 










growing probability in the stochastic region in the decoupled continuous-time mutation- 
selection model can be described using the theoretical formula for the growing probability 
in the Moran two-allele model. 

Figure [3] shows the logarithm of the growing probability, log(Pg), as a function of the 
logarithm of the selective advantage, log(s), for L = 10, p = 0.001 and k = 0.003. In Fig. 
[3l the symbols represent the simulation results for the growing probability in the decoupled 
continuous-time mutation-selection model in a hnite population with a size, N = 10^, for 
a range of measuring parameters, C = 10“^, 10“^, 10“® and 10“®. The parameters in Fig. 
E] belong to the stochastic region, N -C Nc, because the critical population size, Nc, is 
independent of the measuring parameter. Figure [3] shows that the growing probabilities 
approach 1/NC when s* -C 1/NC^ where s* ~ 100s. In Fig. El the theoretical formula 
in Eq. (E]) with N = 10^ and s* = 100s is represented by the short-dotted lines, which 
correspond, bottom to top, to the measuring parameters C = 10“^, lO”"^, 10“® and 10“®. A 
comparison of the symbols and the short-dotted lines shows that the simulation results for 
the growing probability, Pg, can be htted well using the theoretical formula for the growing 
probability, P^, in the Moran two-allele model. 

Figure m shows the logarithm of the growing probability, log(Pg), as a function of the 
logarithm of the selective advantage, log(s), for L = 10, /i = 0.001, k = 0.003, and C = 
10“^. In Fig. m the symbols represent the simulation results for the growing probability 
in the decoupled continuous-time mutation-selection model for a range of population sizes, 
N = 10^, 10® and 10^. The parameters in Fig. E] belong to the stochastic region, N <C Nc, 
as mentioned in the explanation of Fig. [H Figure 0] shows that the growing probabilities 
approach 1/NC when s* -C 1/NC^ where s* ~ 100s. In Fig. 01 the theoretical formula in Eq. 
(1^ with C = 10“^ and s* = 100s is represented by the short-dotted lines, which correspond, 
top to bottom, to the population sizes, N = 10®, 10® and 10^. The simulation results in Figs. 
El El and 0] show that the growing probability, Pg, in the stochastic region in the decoupled 
continuous-time mutation-selection model for various population sizes, selective advantages, 
and measuring parameters can be htted well using the theoretical formula for the growing 
probability, P^, in the Moran two-allele model. 

Figure [5] shows the logarithm of the growing probability, log(Pg), as a function of the 
logarithm of the htness parameter, log(/c), for L = 10, p = 10“®, C = 10“^ and s = 10“^. 
In Fig. El the simulation results for the growing probability in the decoupled continuous- 
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time mutation-selection model are represented by the solid lines with the symbols, which 
correspond, left to right, to the population sizes N = 10®, 10®, 10^® and 10^^. The growing 
probability, Pg, in Fig. [5] represents the saturated growing probability, P*, because C ^ 
1/Ns* for all parameters considered, where s* ~ 100s ~ 10“^. The arrows in Fig. [5]indicate 
the critical htness parameters, kc, left to right, for N = 10®, 10®, 10^® and 10^^, which 
represent the htness parameters satisfying iVAX£°(f = 0) = 1. Figure [5] shows that the 
saturated growing probability in the deterministic region, k <C k^ approaches unity and that 
the saturated growing probability in the stochastic region, k kc, becomes approximately 
the effective selective advantage, regardless of the htness parameter and the population size. 
The simulation results in Figs. El [3], 0] and [5] indicate that the selective advantage ratio, 
r, maintains a similar value for various population sizes, selective advantages, measuring 
parameters, and htness parameters, if the sequence length L is the same. 

Figure [6] shows the logarithm of the growing probability, log(Pg), as a function of the 
logarithm of the selective advantage, log(s), for p = 0.001 and k = 0.002. In Fig. [6l the 
simulation results for the growing probability in the decoupled continuous-time mutation- 
selection model are represented by the symbols, that correspond, top to bottom, to the 
sequence lengths, L = 10, 14 and 18. The computer simulation for an inhnite population 
showed that the critical population sizes satisfying = 0) = 1 were approxi¬ 

mately 3.98 X 10®, 2.48 x 10® and 1.61 x 10^^ for L = 10, 14 and 18, respectively. Hence, the 
population sizes were set to N = 10®, 10® and 10^^ for L = 10, 14 and 18, respectively, to 
belong to the stochastic region, N <C N^. The measuring parameters were set to C = 10“^, 
10“® and 10“® for L = 10, 14 and 18, respectively, to maintain the same product of the pop¬ 
ulation size and the measuring parameter, or NC = 10®, because the growing probability 
depends on iVC', as shown in Eq. ([6]). 

Figure [6] shows that the simulation results for the growing probability for L = 10, 14 
and 18 have a similar dependence on the selective advantage. In Fig. [6l the theoretical 
formula in Eq. ([6]) with NC = 10® and s* = rs is represented by the short-dotted lines, 
which correspond, top to bottom, to the selective advantage ratios, r = 100, 75 and 55. A 
comparison of the symbols and the short-dotted lines showed that the simulation results for 
the growing probability, Pg, could be htted well using the theoretical formula for the grow¬ 
ing probability in the Moran two-allele model with the adjusted selective advantage ratio, 
which decreases with increasing sequence length. The simulation showed that the selective 
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advantage ratio does not depend on the population size, selective advantage, measuring pa¬ 
rameter, and fitness parameter; instead, that the selective advantage ratio decreases with 
increasing sequence length. 


IV. SUMMARY 

The growing probabilities of additional offspring with the beneficial reversal allele for 
growing to a size, iVU, for various population sizes, N, sequence lengths, L, selective advan¬ 
tages, s, fitness parameters, k, and measuring parameters, U, were calculated for a haploid, 
asexual population in the decoupled continuous-time mutation-select ion model using the 
logistic branching process. The growing probability was measured from the initial steady 
state in a sharply-peaked landscape, by switching on the asymmetric sharply-peaked land¬ 
scape with a positive selective advantage of the reversal allele over the optimal allele. The 
growing probability was evaluated as the ratio of the first arrival time to the hnal arrival 
time, where the hrst and the hnal arrival times are dehned as the arrival time of the hrst 
additional offspring and that of a successful offspring with the reversal allele that can grow 
to a size, NC, respectively. 

Computer simulations showed that the growing probability in an asymmetric sharply- 
peaked landscape in the decoupled continuous-time mutation-selection model in the stochas¬ 
tic region satisfying NAX'^{t = 0) <C 1 was inversely proportional to the measuring pa¬ 
rameter when C < 1/Ns*, bent when C ~ 1/Ns* and became saturated when C > 1/Ns*, 
where s* is the effective selective advantage. This suggests that the number of additional 
offspring with the reversal allele, n^, in the stochastic region drifts neutrally until ul reaches 
the size, NC, when ni < 1/s*, and that it grows deterministically when > 1/s*. Com¬ 
puter simulations also showed that the saturated growing probability in the stochastic region 
was approximately the effective selective advantage when C ^ 1/Ns* and s* 1. 

Computer simulations showed that the growing probability in the stochastic region in an 
asymmetric sharply-peaked landscape in the decoupled continuous-time mutation-select ion 
model for various population sizes, selective advantages and measuring parameters could 
be htted well using the theoretical formula for the growing probability in the Moran two- 
allele model. A previous study showed that the growing probability in the stochastic region 
in the coupled discrete-time mutation-selection model could be described using the theo- 
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retical formula for the growing probability in the Wright-Fisher two-allele model[l8|. The 
present study suggests that the growing probability in the stochastic region in an asymmet¬ 
ric sharply-peaked landscape in the decoupled continuous-time mutation-selection model 
could be described using the theoretical formula for the growing probability in the Moran 
two-allele model. 

Computer simulations showed that the selective advantage ratio, which represents the 
ratio of the effective selective advantage to the selective advantage, does not depend on the 
population size, selective advantage, measuring parameter, and htness parameter; instead, 
that the selective advantage ratio decreases with increasing sequence length. The growth 
behavior of additional offspring with the reversal allele in an asymmetric sharply-peaked 
landscape in the decoupled continuous-time mutation-selection model was controlled by the 
effective selective advantage of the reversal allele over the optimal allele despite there being 
many other alleles with lower htness. The present study might improve the understanding 
of the evolutionary process in the decoupled continuous-time mutation-selection model. 
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log(C) 


FIG. 1: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the 
measuring parameter, log(C), for L = 10, = 0.001, k = 0.003 and s = 10“^. The simulation 

results for the growing probability in the decoupled continuous-time mutation-selection model are 
represented by the solid lines with symbols, which correspond, right to left, to the population sizes 
N = 10^ 10®, 10^, 10® and 10®. 
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FIG. 2: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the 
measuring parameter, log(C), for L = 10, /r = 0.001 and k = 0.003. The simulation results for the 
growing probability in a finite population with a size, N = 10^, in the decoupled continuous-time 
mutation-selection model are represented by the symbols, and the theoretical formula in Eq. ([6|) 
with N = 10^ and s* = 100s is represented by the short-dotted lines, which correspond, top to 
bottom, to the selective advantages, s = 10“^, 10“^, 10“® and 10“®. 
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FIG. 3: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the 
selective advantage, log(s), for L = 10, /r = 0.001 and k = 0.003. The simulation results for the 
growing probability in a finite population with a size, N = 10^, in the decoupled continuous-time 
mutation-selection model are represented by the symbols, and the theoretical formula in Eq. Q 
with N = 10^ and s* = 100s is represented by the short-dotted lines, which correspond, bottom to 
top, to the measuring parameters, C = 10“^, 10“^, 10“® and 10~®. 
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FIG. 4: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the selective 
advantage, log(s), for L = 10, /i = 0.001, k = 0.003 and C = lO”'^. The simulation results for the 
growing probability in the decoupled continuous-time mutation-selection model are represented by 
the symbols, and the theoretical formula in Eq. Q with C = 10“^ and s* = 100s is represented 
by the short-dotted lines, which correspond, top to bottom, to the population sizes, N = 10®, 10® 
and 10^. 
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FIG. 5: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the fitness 
parameter, log(A:), for L = 10, /x = 10“^, C = 10“^ and s = 10“^. The simulation results for the 
growing probability in the decoupled continuous-time mutation-selection model are represented by 
the solid lines with symbols, and the critical fitness parameters, kc, are indicated by the arrows, 
which correspond, left to right, to the population sizes N = 10®, 10®, 10^® and 10^^. 
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FIG. 6: Logarithm of the growing probability, log(Pg), as a function of the logarithm of the selective 
advantage, log(s), for fi = 0.001 and k = 0.002. The population sizes are set to = 10^, 10® and 
10^^, and the measuring parameters are set to C = 10“^, 10“® and 10“®, for L = 10, 14 and 18, 
respectively. The simulation results for the growing probability in the decoupled continuous-time 
mutation-selection model are represented by the symbols, which correspond, top to bottom, to the 
sequence lengths, L = 10, 14 and 18. The theoretical formula in Eq. (I6|) with NC = 10® and 
s* = rs is represented by the short-dotted lines, which correspond, top to bottom, to the selective 
advantage ratios, r = 100, 75 and 55. 







